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We study single-flavour quark pairing (“self-pairing”) in colour-superconducting phases of quark 
matter, paying particular attention to the difference between scenarios where all three flavours 
undergo single-flavour pairing, and scenarios where two flavours pair with each other (“2SC” pairing) 
and the remaining flavour self-pairs. We perform our calculations in the mean held approximation 
using a pointlike four-fermion interaction based on single gluon exchange. We confirm the result 
from previous weakly-coupled-QCD calculations, that when all three flavours self-pair the favored 
channel for each is colour-spin-locked (CSL) pseudoisotropic pairing. However, we find that when the 
up and down quarks undergo 2SC pairing, they induce a colour chemical potential that disfavors the 
CSL phase. The strange quarks then self-pair in a “polar” channel that breaks rotational invariance, 
although the CSL phase may survive in a narrow range of densities. 
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I. INTRODUCTION 

Matter at high density and sufficiently low tempera¬ 
ture is expected to form a colour simerconducting con¬ 
densate of quark pairs Q, II i i (for reviews, see 
Ref. Such a phase of matter may exist in 

the cores of compact stars or be created during low 
energy heavy-ion collisions |ll| . At high density, quarks 
will fill up the available energy states to form a Fermi 
surface. Since two quarks in the antisymmetric colour 
antitriplet channel experience an attractive interaction, 
we expect a BGS pairing instability 0 at the Fermi 
surface. The instability is resolved via the opening of 
a energy gap in the quasiparticle spectrum. The quark 
pairing breaks the 5'17(3)coiour gauge symmetry of QCD, 
justifying the name colour superconductor. 

The strong interaction is most attractive between 
quarks in the colour-antisymmetric spin-0 channel, but 
Fermi-Dirac statistics then require flavour-antisymmetric 
pairing, involving two different flavours. Well-studied 
candidates include the 2SG @,0,13 GFL [l^lT^I3| 
phases of quark matter. At asymptotically high density, 
the favored phase of quark matter is the GFL phase. 

However, the cores of compact stars are not asymptot¬ 
ically dense and it becomes necessary to consider real- 
world effects such as a non-zero strange quark mass, elec¬ 
trical and colour charge neutrality and /3-equilibrium. 
Taken together, these constraints act to separate the 
Fermi surfaces of the different quark flavours. This means 
that as density decreases from the asymptotic regime it 
becomes harder and harder to maintain the GFL con¬ 
densate, and we expect transitions to other phases. In 
such a context it is natural to look for single-flavour pair¬ 
ing patterns, and many different ones have been found 


Ellllllliaill. Galculations in weak-coupling QGD 
indicate that for a single isolated flavour the colour-spin- 
locked (CSL) phase has the lowest free energy. How¬ 
ever, in the realistic context of quark matter in neutron 
stars there are three flavours present. In this paper we 
explore the single-flavour phases in such a context, us¬ 
ing a Nambu-Jona-Lasinio model. We find that if the 
up and down quarks undergo two-flavour (2SC) pairing 
then this typically induces a non-zero colour chemical 
potential that disfavors CSL pairing for the remaining 
flavour, which then self-pairs in a different “transverse 
polar” channel instead, which we will simply label as 
“ISC”. 

The I SC / pol ar channel was studied in the NJL con¬ 
text in Ref. |2(Tl| . where we surveyed the pairing patterns 
that factorize into a product of colour, flavour, and Dirac 
(spin) structures (this ansatz excludes the CSL phase, 
where colour and spin are intertwined). We used an NJL 
model with various plausible interactions, and found that 
for a single flavour the most attractive channel was the 
(3^1,35,1,-|-)(C 73 ) channel, where the notation translates 
as (colour antisymmetric 3, flavour symmetric 6, spin 1, 
parity even), with Dirac structure Cya. Pairing between 
quarks of the same flavour in this channel was found to 
have a gap parameter of ~ 1 — 10 MeV. Unlike the CSL 
phase, this condensate breaks rotational symmetry. In 
weak-coupling QCD calculations this is called the “trans¬ 
verse polar” phase |3 • 

In this paper we do not attempt a comprehensive study 
of all the phases that have been suggested for quark mat¬ 
ter. Our purpose is to show that the presence of 2SC 
pairing can affect the single-flavour pairing of the re¬ 
maining flavour, so we restrict our study to polar and 
CSL pairing of individual flavours, and 2SC pairing of 
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FIG. 1: Pictorial representation of simple single-flavour pair¬ 
ing in neutral quark matter. This will be referred to as the 
(ISC)® phase in the text. The requirement of electric neutral¬ 
ity and a nonzero strange quark mass forces the Fermi mo¬ 
menta of the three flavours apart. Two colours of each flavour 
form (3a,3s,1i+)(C'73) Cooper pairs (ISCu, ISCd and ISCs). 
The third colour of each flavour forms (6s,3s,0,-|-)(C7o75) 
pairs. 


FIG. 3: Pictorial representation of CSL pairing in neutral 
quark matter. This will be referred to as the (CSL)® phase in 
the text. It is composed of up, down and strange quark CSL 
condensates, labelled CSLu, CSLd and CSLs respectively. 
The requirement of electric neutrality and a nonzero strange 
quark mass forces the Fermi momenta of the three flavours 
apart. The red, green and blue colours of each flavour pair in 
a colour-antisymmetric channel. 
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FIG. 2: Pictorial representation of simultaneous 2SC and ISC 
pairing in neutral quark matter. This will be referred to as 
the 2SC-|-lSCs phase in the text. This will only occur if the 
condensation energy of the 2SC pairing is strong enough to 
offset the cost of dragging the red and green u and d Fermi 
momenta away from the values dictated by electrical neutral¬ 
ity and the strange quark mass, to a common value. 


FIG. 4: Pictorial representation of 2SC and CSL pairing 
in neutral quark matter. This will be referred to as the 
2SC-|-CSLs phase in the text. Colour neutrality will create 
a small 0{A^/fj.) splitting between the Fermi momenta of the 
blue s quarks and s quarks of the other two colours. If this 
splitting is sufficiently large it will prevent the CSLs conden¬ 
sate from forming. 


the up and down quarks. We do not take into account 
the CFL phase, nor its gapless varia nt g CFL 
crystalline (LOFF) pairing pM and so on. Our 

study is relevant to regions of the phase diagram where 
three flavours of quark are present, but there is either 
no cross-flavour pairing, or only pairing between the up 
and down quarks (“2SC-|-s” in the the nomenclature of 
Ref. m ). Such regions have been found to exist in stud¬ 
ies of the QC D p hase diagram using NJL models. For 
example Ref. finds a 2SC phase with strange quarks 
for stronger diquark coupling {Gd = Gs) at T = 0, over 
the chemical potential range 398 MeV ^ H 'S 412 MeV 
(see their Figs. 2 and 6). Ref. also finds such a phase 
(see their Fig. 7). It is also interesting to note that there 
is some chance our considerations might be relevant to 
the early life of the star, when T > 1 MeV, neutrinos 
are trapped, and the resultant lepton number chemical 
potential can favor a 2SC phase with strange quarks . 
However some studies only find up and down quarks at 
this time [^. Moreover, single flavour phases will only 
be relevant if they have large enough gaps to survive at 
these temperatures. We find that their critical temper¬ 
atures are of order 1 MeV(see Fig. [TJ, but this result is 


sensitive to the cutoff (see end of Section^ so it remains 
an open question whether they can survive during the era 
of neutrino trapping. 

The paper has the following structure. Section ITTIsum- 
marizes the ISC and CSL pairing patterns. Sections cni 
and HU discuss the model and approximations used to 
calculate the free energy of the colour superconducting 
phases that are of interest here. In all cases, electrical 
and colour neutrality constraints are imposed. Section Ivl 
shows the behaviour of the ISC and CSL gap parameters 
as functions of ^ and T. The gaps are then used in a free 
energy calculation, the results of which are presented in 
Section EH as a sequence of (/r,r) phase diagrams for 
neutral three flavour quark matter for different values of 
Mg. Section Em presents conclusions and directions for 
further research. 
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II. ISC AND CSL PAIRING 

The forms of zero-temperature pairing that will be 
compared in this paper are shown in Figs. ^ 121 121 
m In the “(1SC)3” phase (Fig. UJ, the Fermi mo¬ 
menta of the three flavours are so different that only 
single-flavour pairing is possible. Two colours of each 
flavour form (3^,35,1,-|-)(C'73) Cooper pairs (ISCu, 
ISCd and ISCs). The third colour of each flavour 
forms (6 s,3s,0,-|-)(C'7o 75) pairs. In the 2SC-|-ISCs phase 
(Fig. El, the up and down quarks are able to pair with 
each other in the 2SC channel (which only involves two 
colours), but the strange quarks undergo single-flavour 
pairing. Again, the third colour of each flavour forms 
(6s,3s,0,-f)(C7o75) pairs. 

Figures 01 and 0 are the equivalents of Figures Hand 01 
where the ISC pairing of individual flavours is replaced 
by CSL pairing. 

We now give the explicit form of the pairing patterns 
mentioned so far. The basic two-colour single-flavour 
spin-1 ISC pairing (the “polar phase” in Ref. 0 ) , de¬ 
scribed here as (3 a,3s, 1,-|-)(C73), is 

= Aisc(C73)“'’(A2)a/3<5*^' . (1) 

The colour indices are a, /3, and A 2 is an antisymmet¬ 
ric Gell-Mann matrix in colour space. The Dirac indices 
are a, b, and C = 17072 is the usual charge conjugation 
matrix. The flavour indices i,j are included for com¬ 
pleteness: in this single-flavour channel i = j = \. This 
condensate breaks rotational invariance, as the explicit 
occurrence of the 73 matrix indicates. Only two colours 
(typically taken as red and green) of quark are involved 
in the pairing meaning that the remaining (blue) quark is 
required to find an alternative channel in which to pair. 
The S'f7(3)coiour symmetry is broken down to an SU(2) 
subgroup. 

The single-colour single-flavour ( 65 , 3 s, 0 ,-|-)(C' 7 o 75 ) 
pairing is 

= A05(C'7075)“'^a/35*'' . (2) 

This channel is attractive in an NJL model with gluon- 
exchange-type interaction, but the pairing disappears as 
the quark masses go to zero [^. This means it may be 
somewhat suppressed for up and down quarks, depending 
on how much of their vacuum constituent mass survives 
at the relevant densities [n m 03 . Our earlier NJL 
study estimated the pairing gap in the Cjoj5 channel 
to be of the order 1-100 eV. This is much smaller than 
the typical gaps in the 2SC, ISC, and CSL channels, so 
when we perform free energy calculations we will neglect 
the ( 6 s, 35 , 0 ,-|-)(C' 7 o 75 ) pairing, treating those quarks as 
unpaired. 

CSL pairing is 

=AcsL{CjAr^\if,S^\ (3) 

where we sum over A = 1,2,3, 7 ^ = ( 71 , 72 , 73 ) and 
Xa = (A 7 , A 5 , A 2 ) which are the antisymmetric Gell-Mann 


matrices in colour space. The symmetry breaking asso¬ 
ciated with the CSL phase is 

5C(3)colour X SOi3)j X U{1)b ^ 50(3)eolour+J X Z 2 , 

(4) 

where S'f7(3)coiour is gauged, SO{3)j is the rotation 
group and U{1)b is baryon number. The unbroken rota¬ 
tional symmetry S'0(3)coiour+j corresponds to a locking 
between the colour and spin degrees of freedom. This can 
be thought of as associating a direction in configuration 
space with a linear combination of directions in colour 
space, i.e. red and green quarks pair in the z-direction, 
red and blue quarks in the y-direction and green and 
blue quarks in the x-direction. Therefore, unlike the po¬ 
lar phase, CSL pairing remains isotropic. The locking 
is similar to that which is found in the CFL phase be¬ 
tween colour and flavour. The CSL phase is the colour 
superconducting equivalent of the R-phase of superfluid 
helium-3 Q. 

In general, one would expect that the CSL conden¬ 
sate could involve a linear combination of Cjt and Caoi 
with a colour structure that is correlated with the spatial 
direction. However, Table 1 in Ref. shows that the 
colour antisymmetric Cctos channel gives no contribution 
when using the mean-field approximation within an NJL 
model with a full gluon interaction vertex. This is pre¬ 
cisely the interaction that we use, so only the Cji term 
is included in the CSL ansatz of Equation m- 

We can calculate the binding energy of the CSL chan¬ 
nel, following Ref. |^, and compare it with the results 
for the polar ISC channel, taken from the same refer¬ 
ence. The result is shown in Table ^ where for a given 
condensate of gap parameter A the binding energy via 
interaction I is given hy E = — S/A^, so larger values 
of S imply stronger binding. / can be “inst” (instan- 
ton interaction), “mag” (magnetic gluon interaction), or 
“elec” (electric gluon interaction). The table shows that 
the CSL condensate is more tightly bound than the polar 
condensate, leading one to expect that, other things be¬ 
ing equal, CSL pairing will be favored over polar pairing. 


III. OUTLINE OF CALCULATION 

The aim of this paper is to map out the (^,T) 
phase diagram for the 2SC, unpaired, (ISC)^, (CSL)^, 
2SC-|-lSCs, and 2SC-|-CSLs phases of neutral three 
flavour quark matter at densities relevant for compact 
stars (^ ~ 500 MeV), for different values of Mg. The 
phase diagram is obtained by calculating the free energy 
of each phase, imposing neutrality conditions and equili¬ 
bration under the weak interactions, and minimizing the 
free energy with respect to the pairing strengths Ap in 
the various phases. In section ITVI we will discuss the NJL 
model within which the free energy is calculated. In this 
section we assume that the free energy is known. 

To obtain neutral matter we introduce chemical po¬ 
tentials coupled to the gauged charges. (In full QCD 
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phase 

colour 

flavour 

j 

parity 

Dirac 

BCS 

enhance¬ 

ment 

‘S'inst 

5*elec 

H” 5*mag 

‘5'mag 

CSL 

3a 

3s 

Is 

-f 

Cyi LR 

0(1) 

0 

-t96 

-k48 

“Polar” (ISC) 

3a 

3s 

Is 

-f 

C 73 LR 

0(1) 

0 

-t32 

-fl6 


TABLE I: Binding strengths of CSL (1st row) and ISC (2nd row) channels in an NJL model in the mean-field approximation. 
The ISC results are taken from Table 1 in We see that CSL is more strongly bound. 


the time-components of the gauge fields play this role 
1 ^ . and there are no gluonic contributions to the 
charge density [13 ■) The electron chemical potential fie 
couples to negative electric charge. To impose colour 
neutrality it is sufficient to couple chemical potentials ^3 
and fis to the Cartan sub-algebra of the colour gauge 
group, generated hy J\= diag(l/2, —1/2, 0) and Tg = 
diag(l/3,l/3,-2/3) (H. 

The charge densities are related to the chemical poten¬ 
tials via the derivatives of the free energy, so the neutral¬ 
ity constraints are 


718 = - 


713 = 


77Q = 


dn 

dfis 

dn 

9/i3 

dn 

dfle 


0 , 

0 , 

0 . 


(5) 


The weak interaction violates flavour conservation but 
conserves baryon number, allowing the chemical poten¬ 
tials of the individual quark species to be written in terms 
of lie, fi 3 , Ms and the quark number chemical potential 

M, m 

i-^ru — 3 M 85 

IJ^gu = //- |/Xe - ^M3 + 

2 2 

f^bu — 3 /^^ 3 M 85 

Mrd = M + \^J‘e + + |M8, 

^gd ^ 3/^6 2^2 3M8; (0) 

12 ^ ^ 
f^bd — M “I” 3 /^^ 3 M 85 

I 1 I 1 I 1 

lirs = M+3Me+2M3+3M8-^-^, 

I 1 111 Af/ M/ 

figs = M+3M«- 2/^3+3M8-^-^, 

_ I 1 2 

Mbs — li+ ^ fie ^flS 2fj. 8/^3- 


The strange quark mass is included to leading order as a 
term in the strange quark chemical potential IHiii. 

The neutrality constraints will always be implemented 
locally, i.e. no mixed phases of q uark matter and nuclear 
matter will be discussed [^ [^ . 

For each condensate p, the gap equation is found an¬ 
alytically by differentiating the free energy with respect 
to the gap parameter of that condensate. 


~dA;r 


(7) 


Finding the free energy of a phase therefore corresponds 
to finding a stationary point of n in the parameter space 
of chemical potentials and gap parameters. If we work in 
the neutral subspace, fixing Me, M 3 j Ms as functions of the 
Ap, then the relevant stationary point will always be a 
minimum with respect to the gap parameters. 

We can now anticipate some simplifications that will 
streamline our calculation. 

Firstly, as is clear from Figs, mill El El we will always be 
concerned with colour superconducting condensates that 
are at least invariant under the SU(2) colour subgroup 
that rotates red and green quarks into each other. This 
means that M 3 will always be zero. 

Secondly, it is never necessary to solve the full set of 
coupled gap and neutrality equations. The single-flavour 
pairing is so weak that its back-reaction on the chemical 
potentials via the neutrality condition can be neglected. 
Thus for the 2SC-|-lSCs and 2SC-|-CSLs phases one can 
first treat the strange quarks as free, determine Me and 
Ms by solving for 2SC pairing, and then use those values 
in the gap equations for the strange quark pairing. Sim¬ 
ilarly, for the (ISC)^ and (CSL)^ phases one can use the 
values of Me and Ms that would ensure neutrality in to¬ 
tally unpaired quark matter (in the case of Ms this value 
is zero). 

Actually, the only phase where the exact value of Ms is 
important is the 2SC-|-CSLs phase, where the nonzero Ms 
puts a stress on the CSLs pairing. Following the standard 
argument [^, we expect that CSL pairing can occur as 
long as the cost of breaking a CSL pair is greater than 
the gain in energy from converting a blue quark to a 
red/green one, i.e., 


|Ms|<2Acsl- (8) 

The phase diagrams to be presented in section EH will 
show this behaviour. 


IV. MODEL AND APPROXIMATIONS 

The NJL-type Lagrangian that will be used in the fol¬ 
lowing calculations is, 

^NJL = -Ckin + -C int ( 9 ) 

= 7/(7^- + Ant, ( 10 ) 
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where jji = /i 7 o is a diagonal matrix in colour-flavour 
space, with entries given by 0 . The up and down quarks 
are treated as massless, and the strange quark mass is 
only taken into account to leading order, as an effective 
chemical potential for strangeness. The interaction we 
use is abstracted from single-gluon exchange by replac¬ 
ing the gluon propagator with a four-fermion coupling 
constant, 

Ant = -^G(^r“V’)(^r^V'), (11) 

where 


r“ = (12) 

where the A“’s are the standard Gell-Mann matrices (a = 
1.. .8), giving the interaction the quantum numbers of 
single gluon exchange. 

In order to study diquark pairing and colour super¬ 
conductivity, we Fierz transform the interaction into the 
form 

= (13) 

p 


where p labels the different diquark pairing channels, T®"^ 
specifies the form of the diquark pairing in channel p and 
f 99 is defined by 

rr=7o(f«)t7o^- (14) 

Note that we sum up the terms for all the different 
pairing channels, even though no single Fierzing of the 
original interaction would give the sum of all of them. 
This is legitimate because there are no cross-terms be¬ 
tween different channels in the binding energy. 

With the interaction expressed as in Eq. Ha it is 
straightforward to do a mean field calculation of the free 
energy is now in a suitable form to allow the free energy 
to be calculated. This is done by performing a bosonisa- 
tion of the Fierz transformed La grangia n via a Hubbard- 
Stratonovich transformation [42.l4dll44Lliq . The general 
expression for the thermodynamic potential of a colour 
superconducting condensate is 


n = -T 


^ J (27r)3 2 \ T ) Aa'G ’ 

(15) 

where S ^{iujn,p) is the full inverse quark propagator. 
ujn = (2n — l)7rT are the Matsubara frequencies and the 
trace is over colour, flavour, and spinor indices. The 
Matsubara summation is performed using the identity 


r^ln(' ^' +^a(pT) \ ^ |,^(p-)|+2Tln(l+e-l^“(^^l/^). 

n ^ ^ 

(16) 

The functions ea(p) are the dispersion relations for the 
fermionic quasiparticles. They are not explicitly T- 
dependent, but they depend upon the gap and quark 
chemical potentials which are T-dependent. 


For convenience, we use neutral unpaired quark mat¬ 
ter as the zero of free energy, so we subtract it from the 
value computed above. This gives a physically meaning¬ 
ful quantity without ultraviolet divergences: 


n = - 


-f 



eo(^| + 2rin(l-I-e 


(17) 

where the electron contribution to the free energy has 
been included. The electron chemical potential in neutral 
unpaired quark matter is 


. (18) 

The quasiquark dispersion relations ea(p) are the values 
of the energy at which the propagator diverges, i.e. 

det5'“^(ea(^,^ = 0. (19) 


The calculation of then depends upon the calculation 
of the determinant of the matrix S~^. This matrix can 
be block-diagonalized in the colour-flavour space, with 
one block for each independent pairing channel. For 
each colour and flavour there are eight dispersion re¬ 
lations (four Dirac branches including left-handed and 
right-handed, particle and antiparticle, doubled by the 
Nambu-Gorkov formalism). In calculating the free en¬ 
ergy of quark matter with three flavours and three colours 
the index a in P|) will therefore vary from 1 to 72. 


A. Dispersion relations 

For each unpaired colour and flavour of quark, there 
two branches of the dispersion relation (particle and an¬ 
tiparticle) each with multiplicity 4: 

^Le = P ± (20) 

For each condensate we now give the dispersion rela¬ 
tion of the quasiquarks (with their multiplicity) and the 
binding energy parameter a'. 


1. 2SC 


For 2SG pairing of the red and green quarks of up and 
down flavour, the dispersion relations are 


i?2SC = \/(P ± p)^ -I- A2 ± Sp 


( 21 ) 


each with multiplicity 4 (spin and Nambu-Gorkov), 
where 


M — Q {Pru T Pgd') — {Pru 4” pgd') 

Pe , P8 

= T 

1 6 3 1 ^ 

Sp = l^^pgd ~ Pru) = '^ip'rd ~ pgu) = 


( 22 ) 
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The binding energy parameter is 


V. RESULTS 


(23) 


If A 2 SC < 2SC condensate is in the gapless 

2SC state |Sli3- 


2. ISC 

For ISC pairing of red and green quarks of some flavour 
/ (assumed to have common chemical potential /l), the 
dispersion relations are 

^isc = p'^ + + A'^ ± 2^p?{p^ + m?) + , 

(24) 

each with multiplicity 4. The binding energy parameter 
is 


(25) 


3. CSL 


For CSL pairing of three colours of a single massless 
flavour, in the case ps = 0 , we can obtain a simple closed- 
form expression for the determinant of the inverse prop¬ 
agator, 

det S~~^{po, p) = 

{A^ - {p-poff 

X {A’^ + - {p + pqYY 

X Ua'^ + p^ + (p2 - p2)2 _ + pI) 

X 2 (26) 

-fA^(5^^ -f {3p - 5po){p + po))j 

X (aA^ + + {p'^ - pI)'^- 2p,'^{p‘^+ pI) 

+A'^{5p'^+ {3p + 5po){p-po))^ , 


which has dimension 24, as expected. From this the dis¬ 
persion relations of the 24 quasiparticles can be obtained. 
Note that there are 4 gapless modes (first line of the de¬ 
terminant), and the remaining 20 are gapped. When we 
include a nonzero quark mass all the quark modes are 
gapped, which is different from the behaviour in weak- 
coupling QCD 1^ . For a discussion of this point see 
Refs. 1^ 1^. In the general case ps 7 ^ 0,m > 0, the 
determinant can be readily obtained using a symbolic 
mathematics program such as Mathematics but is too 
complicated to display here. 

The binding energy parameter is 



The following sections present the results of solving 
the gap and neutrality equations and calculating the 
free energies for the (ISC)^, (CSL)^, 2SC-|-lSCs, and 
2SC-|-CSLs condensates. The momentum cutoff was 
fixed at A = 600 MeV, with the value of the coupling G 
calibrated to give A 2 sc{p = 500, Mg = 0) = 31.1 MeV, 
as in Ref. |^. This corresponds to AcFL(Ai = 500, Mg = 
0) = 25 MeV. We also performed some calculations at 
A = 800 MeV with the same calibration, to check for 
cutoff sensitivity. 

If single-flavour pairing were ignored, the phase dia¬ 
gram would contain a transition from unpaired quark 
matter to 2SC quark matter at p = M‘^/(2p). In princi¬ 
ple this critical value also depends on temperature, but 
at the temperatures of interest to us (T < 1 MeV) this 
can be neglected. Thinking of this as a background, 
we can then include the possibility of ISC and CSL 
single-flavour pairing, and we expect to find that in the 
low-temperature parts of the unpaired region there are 
regions of (ISC)^ or (CSL)^ phases, and in the low- 
temperature parts of the 2SC region there are regions 
of 2SC-|-CSLs or 2SC-|-lSCs. This is indeed what we 
ultimately find. 

As mentioned in section CYI we can calculate pe and 
/ig in the 2SC/unpaired background, before including the 
effects of single flavour pairing. We show these results in 
the next subsection. Then we go on to solve the gap 
equations for the various forms of single flavour pairing, 
and compare their free energies. 


A. Chemical potentials with no single flavour 
pairing 

Figure 0 shows (for four different strange quark 
masses) the values of pe and /ig in neutral unpaired/2SC 
matter at T w 0 over a range of quark chemical potential 
450 MeV < p < 550 MeV. This required simultane¬ 
ously solving the 2SC gap equation and the two neutral¬ 
ity equations, and choosing the phase with the lowest free 
energy at each value of p. 

At /i « Mg/(2/i) we see the transition from unpaired 
to 2SC, at which pe jumps from about Mg/{Ap) to about 
Mg/{2p) [^. In the unpaired phase, /ig is zero because 
nothing picks out any direction in colour space. In the 
2SC phase, however, the red and green light quarks pair, 
and this means that a nonzero /ig is required to maintain 
colour neutrality. At moderate densities /ig is negative 
because the number of red and green quarks is enhanced 
by the 2SC pairing (there is more phase space above the 
Fermi surface than below it), and a negative value of 
/ig restores the colour balance by favoring blue quarks 
(Tg = diag(l, 1, —2)). As the quark chemical potential 
approaches the cutoff, however, the situation is reversed: 
degrees of freedom above the cutoff are excluded, so the 
condensate contains more red and green holes below the 
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FIG. 5: For 2SC/unpaired quark matter, the variation of ns 
with quark chemical potential at T « 0 and for various values 
of Ms- We see the transition from unpaired {ns = 0) at low 
n to 2SC at high n- 

Fermi surface than red and green quarks above it. Thus 
ns turns around and starts to rise again as n approaches 

A. Depending on the values of the strange quark mass 
and the cutoff, ns may cross through zero and become 
positive at some specific density, as seen in Fig. El This 
is a very interesting phenomenon, because it leads to 
an “island” of 2SC+CSLs pairing in the phase diagram 
(Fig. El, but we emphasize that it is sensitive to the value 
of the cutoff, and if we work at A = 800 MeV then we 
find that ns only becomes positive at /r > A. 


B. Gap parameter plots 

We now proceed to solve the gap equations for the 
single-flavour-pairing phases. Figure El shows how the 
single-flavour pairing gap parameters vary over the quark 
chemical potential range 450 MeV < n < 550 MeV. In 


all cases, the gaps for the different flavours follow the 
same pattern. 

For the purely single-flavour phases ((ISC)^ and 
(CSL)^) there are three gap parameters, for up, down, 
and strange quarks (red and green only, see Figs. Illdll . 
These are shown by the various broken lines in Fig. El 
They exhibit the hierarchy, 

> A“ > AV (28) 

This follows from the electrical-neutrality-induced or¬ 
dering of Fermi momenta in unpaired quark matter 
Pp>p'P> Pp and the fact that quarks with larger Fermi 
momentum have more phase space at their Fermi surface, 
and hence stronger pairing. We also see that ISC pair¬ 
ing gives a larger gap than CSL pairing, but it will turn 
out that CSL has a larger condensate binding energy (see 
Table 0 and so, other thing being equal, we expect CSL 
to be favoured over ISC. This is confirmed by the free 
energy calculation and can be seen by the presence of the 
CSL phase in the phase diagrams of section IVTI 

In the 2SC-|-CSLs and 2SC-|-lSCs phases only the 
strange quarks undergo single-flavour pairing, so there is 
only one gap parameter. A®. This is shown by the solid 
lines in Fig. El- In the top left panel of the hgure, the 
solid line (2SC-|-lSCs) ends at /r « 510 MeV, because the 
2SC phase ceases to be favored below that point. Notice 
that in general the value of A® in the presence of 2SC 
pairing of the up and down quarks is slightly higher than 
its value when all three flavours undergo single-flavour 
pairing. This is because 2SC pairing pulls the up and 
down Fermi surfaces together, which drags the strange 
Fermi momentum up a little in order to maintain elec¬ 
trical neutrality. This increases the phase space at the 
strange quark Fermi surface, boosting the gap by a small 
amount. 

The gap parameters for the CSLs condensate in the 
2SC-|-CSLs phase show a non-monotonic behaviour (solid 
lines in lower two panels of Fig. El- That is because CSL 
pairing is very sensitive to the value of the ns chemical 
potential, which tries to split apart the Fermi momentum 
of the blue quarks from those of the red and green quarks, 
and the background 2SC pairing induces a non-zero ns 
to maintain colour neutrality. CSL pairing is disfavored 
when its gap is less than a critical value of order ns- For 
the cutoff that we used, A = 600 MeV, ns goes through 
zero at ^ ~ 530 MeV (see Sect. EZl, so the CSL gap 
is boosted there. For a different cutoff, this would have 
occurred at a different density. 

C. Non-BCS temperature dependence 

The temperature dependence of the single-flavour gap 
parameters is shown in Fig. 0 We see that the gap 
parameters drop with increasing temperature, and that 
each vanishes at a second-order phase transition at some 
critical temperature Tent- The critical temperatures de¬ 
viate slightly from the BCS prediction Tcrit — 0.57 A(T = 



















(1SC)3 and 2SC+1SCs 
M=250 MeV M=175MeV 

s s 



(CSL)3 and 2SC+CSLs 
M=250 MeV M=175MeV 

s s 



FIG. 6: Variation of the ISC (upper) and CSL (lower) gaps 
with quark chemical potential at T = 0.03 MeV and for Ms = 
250,175 MeV. Solid lines are for the strange quarks, in the 
case where up and down quarks are undergoing 2SC pairing. 
The various broken lines are for the three flavours of quark in 
the case where all of them undergo single-flavour pairing. 


0 ): 


ISC : 

crit 

a 0.51A{sc(T = 

= 0), 

(29) 

CSL : 

Tf p 

crit 

a 0.82A^SL(r = 

= 0). 

(30) 


Such deviation was also found in weak-coupling-QCD 
studies of single-flavour pairing |^, and non-BCS crit¬ 
ical temperatures have also been seen in gapless phases 
and crystalline phases [^ . 

It is interesting to note that because each flavour has 
a different critical temperature for its self-pairing, we ex¬ 
pect a series of phase transitions as the system cools 
(imagining these phases to be in the core of a neutron 
star, for example), with single-flavour pairing appearing 


(1SC)^ and 2SC+1SCs 
M =250 MeV M =175 MeV 

s s 



(CSL)3 and 2SC+CSLs 
M =250 MeV M =175 MeV 

s s 



FIG. 7: Variation of the ISC (upper) and CSL (lower) gaps 
with temperature at /r = 540 MeV and for Ms = 250,175 
MeV. Solid lines are for the strange quarks, in the case where 
up and down quarks are undergoing 2SC pairing. The various 
broken lines are for the three flavours of quark in the case 
where all of them potentially undergo single-flavour pairing. 


first for the down quarks, then for the up quarks, and 
Anally for the strange quarks. 

In the absence of any 2SC pairing, we also see that 
the critical temperature for ISC pairing of any given 
flavour is approximately equal to the critical tempera¬ 
ture for CSL pairing of the same flavour, even though 
the ISC phase has a larger gap. This means that if we 
heat up a (CSL)^ phase, we do not expect it to turn 
into (ISC)^. However, in the presence of 2SC pairing, 
we find that the critical temperature for 2SC-|-CSLs is 
smaller than that for 2SC-|-lSCs, so we expect to And 
2SC-l-CSLs—>2SC-l-lSCs—>2SC as we heat up the system, 
and this is seen in the phase diagrams of section IVll 

We repeated our calculations for a higher cutoff A = 







































9 


M =250 MeV M =225 MeV 

S S 





FIG. 8: Phase diagrams for one- and two- flavour pairing in nentral three flavour quark matter with values of Ms between 175 
and 250 MeV, calculated in an NJ model with cutoff A = 600 MeV. For the 2SC+lSCs and 2SC-|-CSLs phases see Figs. Inland 
0 The phase labelled “CSL” involves CSL pairing of some or all of the u,d, and s quarks. 


800 MeV, calibrating the coupling to give the same 2SC 
gap, and we found that the single-flavour pairing gaps 
dropped by a factor of ^ 8, but the relations were 
still valid. 


VI. DISCUSSION OF PHASE DIAGRAMS 

Section nvi describes how we calculate the free energy 
for each phase as a function of (/r, T, Ms). At each value 
of (ir,T) the favored phase is the one with the lowest free 
energy. Fig. |S1 presents this information as a series of 
phase diagrams. 

The crude structure of the diagrams involves the ex¬ 
pected transition from unpaired quark matter at low den¬ 
sity to two-flavour-paired quark matter (2SC) at high 
density, at a critical chemical potential /icrit ~ 

In the ^-T diagram this transition occurs along a line that 
comes down very steeply at /i = ^crit, because we are in¬ 


terested in low temperatures T < 1 MeV, so T ^ ^ 2 SC- 
As the strange quark mass drops, this line moves to the 
left, and 2SC occupies more and more of the phase space. 

In both the unpaired and 2SC regions of phase space, 
we find that, at sufficiently low temperatures, single¬ 
flavour pairing occurs. Below ^crit at zero temperature, 
the quark matter is in the (CSL)^ phase in which all 
three flavours undergo CSL pairing, as depicted in Fig.|3 
Since the three flavours have different Fermi momenta, 
they have different CSL gaps, and different critical tem¬ 
peratures, so the region we have labelled in the phase 
diagrams as “CSL” actually contains three bands, one 
where just down quarks undergo CSL pairing, then below 
that one where down and up quarks undergo CSL pair¬ 
ing, and finally at the lowest temperatures one where all 
three flavours undergo CSL pairing. We do not expect 
any ISC pairing in this region because the ISC phases 
do not have higher critical temperatures than the CSL 
phases (Fig. [TJ. 
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Above /icrit, the phase structure is more complicated. 
A colour chemical potential /ig is generated in order to 
keep the 2SC phase colour neutral, and this acts to split 
the blue quark Fermi momentum apart from the red and 
green quark Fermi momenta. The value of /ig is typi¬ 
cally of order a few MeV, which is easily large enough 
to disrupt any attempt at CSL pairing, which itself only 
has a gap on the order of one MeV. As a result, the 
favored phase at low temperature for /i > /Xcrit is typi¬ 
cally the 2SC-|-lSCs phase, in which the red and green 
strange quarks form a spin-1 condensate (Fig. [21). At a 
particular value of fi, however, //g may pass through zero, 
and in that region of phase space, where |^g| < 2Acsl, 
the strange quarks of all three colours are able to pair 
in a 2SC-l-CSLs phase (Fig. which produces an “is¬ 
land” of 2SC-|-CSLs in the 2SC-|-lSCs band. For our 
particular choice of NJL model and cutoff, this occurs 
at ^ « 530 MeV for Ms = 175, 200, 225 MeV. As we 
emphasized in section Sect. jvH however, the position of 
such an island is sensitive to details of the cutoff, and is 
not a robust prediction of the model. 


VII. CONCLUSIONS 

In this paper we have shown that single flavour pair¬ 
ing is strongly affected by what is happening to other 
flavours. In particular, two-flavour pairing of the light 
flavours induces large enough colour chemical potentials 
to change the nature of the favored single-flavour pair¬ 
ing pattern for the strange quarks. Our results, shown 
in Fig.jSl were obtained using a NJL model in the mean 
field approximation with a pointlike four-fermion inter¬ 
action whose index structure is based on single gluon ex¬ 
change. We treated the up and down quarks as massless, 
and treated the strange quark mass as a fixed parameter. 
Some of the interesting features, such as the 2SC-|-CSLs 
“island”, are cutoff-dependent, but we found that the 
basic structure of the diagram is a robust prediction of 
the model. At high temperatures (T ^ 1 MeV), only 
flavour-antisymmetric pairing such as 2SC can survive. 
At sufficiently low temperatures we find either CSL pair¬ 
ing of some or all flavours, or 2SC-|-lSCs pairing. The 
critical temperature for the formation of these conden¬ 
sates is not reliably predicted by our model, but appears 
to be of order 1 MeV. 

Our analysis was focused on the 2SC/1SC/CSL com¬ 
parison. We ignored other single-flavour pairing patterns 
such as those with Dirac structure C/jq^ and Cyoys, since 
for an NJL model with our interaction they have much 
smaller gap parameters, and hence either make a neg¬ 
ligible contribution to the free en ergy , or are disfavored 
relative to the ISC (Cya) phase |^. (It is interesting 
to note that weak-coupling QCD calculations find that 
a transverse planar phase is the strongest competitor to 
CSL [23; this phase has not yet been studied in the NJL 
context.) More importantly, we also explicitly ignored 
other possible competing phases such as the CFL phase 


and crystalline (LOFF) pairing, and we also ignored chi¬ 
ral condensates that would give larger (density depen¬ 
dent) quark masses. One natural direction in which to 
extend this work is to perform a simultaneous comparison 
of all known phases, along the lines of Refs. al¬ 

though a definitive calculation is not yet possible because 
we still do not know what tr ue g round state underlies the 
unstable gapless phases 1471IH^. l5^ . 

Another topic for future work is the possible ramifi¬ 
cations of our results for the phenomenology of quark 
matter in neutron stars. The phase diagrams of Fig. jS] 
show that a core of quark matter inside a neutron star 
could have a complicated radial structure of different 
phases, each with its own special properties. Moreover, 
this structure would change over time as the star’s inte¬ 
rior cooled. In our calculation the single-flavour phases 
had gaps of order 1 MeV, so the phase structure would 
settle down less than a year after the supernova, but 
it should be noted that we completely ignored the very 
weak single-flavour single-colour pairing that accompa¬ 
nies most of these phases, whose gap parameter could be 
as small as 1 eV [^ , leading to the emergence of further 
phase structure after millions of years. 

For neutron star phenomenology, we are interested in 
the transport properties of the different phases, in par¬ 
ticular their permeability to magnetic fields, as well as 
thermal properties, neutrino emissivity and mean free 
path, viscosity, and so on. These are sensitive to the low- 
energy degrees of freedom, such as Goldstone bosons and 
ungapped quark modes. We postpone a detailed study 
and now sketch some of the obvious features. First, some 
general points: a massive flavour with CSL pairing is 
completely and isotropically gapped with an energy gap 
of approximately mA//i for m, A ^ /x, so some of the 
quasiquarks become gapless in the massless limit (see dis¬ 
cussion after Eq. ra i. A massive flavor with ISC pairing 
has a direction-dependent gap, and in the massless limit 
becomes ungapped at two points on the Fermi surface. 
However we must remember that, like 2SC pairing, ISC 
pairing only involves two of the three colours, so there 
may be other gapless modes involving the third colour. 
Focussing on the phases that we have discussed here, the 
most completely gapped phase is (CSL)^, in which all 
quark modes are gapped, although if the up and down 
flavours have low constituent masses there will be modes 
with rather small gaps. (From Ref. j^. Fig. 4, we ex¬ 
pect niu, md ~ 20 MeV at /x ~ 400 MeV, so assuming a 
CSL gap parameter of 1 MeV we find an energy gap of 
about 50 keV.) Next comes the 2SC-|-CSLs phase, which 
has blue up and down quarks that are effectively gapless 
(they may have an eV-scale gap due to pairing in some 
very weak channel such as Cjols, which is suppressed 
for light quarks). Finally, the 2SC-|-lSCs phase has in 
addition the blue strange quarks (which can generate a 
bigger Cjqj^ self-pairing gap because the strange quark 
is heavier). The 2SC-l-lSCs phase is also special because 
rotational invariance is broken by the ISC pairing, and 
for the red and green strange quarks the gap varies over 
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the Fermi surface, but is always of order [^. The 

full cooling phenomenology of these phases has not yet 
been worked out, although such questions have already 
been discussed for isolated single flavours 
two-flavour quark matter . 

The magnetic properties of these phases are also ex¬ 
pected to show some variety. In unpaired quark matter 
there is no Meissner effect, so magnetic flux is not ex¬ 
pelled or confined to flux tubes. In an idealized 2SC 
phase the four participating quark species are uniformly 
gapped, but the remainder are gapless. The colour SU (3) 
and electromagnetic C/(l) gauge symmetries are broken 
down to a colour SU (2) (giving masses to five of the glu¬ 
ons) and a rotated electromagnetic U{1)q. In the (CSL)^ 
phase all the gauge bosons have Meissner masses. The 
gluon masses are of order gAcsh, and the photon mass 
is of order eAcsL liiEl. The magnetic properties of 
the 2SC-|-CSLs phase have not yet been calculated: we 
expect mixing of the photon and the eighth gluon, with 
Meissner masses that depend on g as well as e. In the 
2SC-|-lSCs phase, under the U{1)q which survived the 
formation of the 2SC condensate, the red, green, and blue 
strange quarks have charges —1,0,0 respectively, so the 


additional ISCs pairing breaks the U(1)q symmetry, giv¬ 
ing the Q photon a mass of order eA^g^. The gluons have 
much larger masses of order gA 2 sc- It remains to be seen 
whether this phase is a type-I or type-II superconductor 
for the Q photon. 

We conclude that although single-flavour pairing 
makes a very small contribution to the free energy of 
quark matter, it contributes interesting extra structure 
to the low-temperature part of the quark matter phase 
diagram. Single-flavour pairing determines the symme¬ 
tries and low-energy-excitations in this region, making 
it directly relevant for the study of signatures of quark 
matter in neutron stars. 
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